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Abstract. In this note we give a new method for getting a series of approxi- 
mations for the extinction probability of the one-dimensional contact process by 
using the Grobner basis. 

1 Introduction 

Let X = {0, l} zd denote a configuration space, where Z d is the <i-dimensional 
integer lattices. The contact process {r] t : t > 0} is an X- valued continuous- 
time Markov process. The model was introduced by Harris in 1974 pQ and 
is considered as a simple model for the spread of a disease with the infection 
rate A. In this setting, an individual at x G Z d for a configuration rj G X is 
infected if r](x) = 1 and healthy if r](x) = 0. The formal generator is given 
by 

x&L d 

where rf G X is defined by i] x (y) = rj(y) (y ^ x), and t] x (x) = l — r](x). Here 
for each x G Z d and r\ G X, the transition rate is 

c(x,r}) = (1 -T)(x)) x A ^2 V(y)+V(x), 

y:\y-x\=l 
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with |x| = \x\ | + • • • + \xa\ ■ In particular, the one- dimensional contact process 
is 



001 - 


-> on 


at rate 


A. 


100 - 


-> 110 


at rate 


A, 


101 - 


-> 111 


at rate 


2A, 


1 - 


-> 


at rate 


1. 



Let Y = {A C U 1 : \A\ < oo}, where \A\ is the number of elements in A. 
Let £^(c Z d ) denote the state at time t of the contact process with ^ = A. 
There is a one-to-one correspondence between £^(c Z d ) and r) t E X such 
that x G £^ if and only if rjt(x) = 1. For any A 6 7, we define the extinction 
probability of A by lim^oo P(£^ = 0). Define v\(A) = v x \r\ : ^(x) = 
for any x E A}, where v x is an invariant measure of the process starting 
from a configuration: r\{x) — 1 (x G Z d ) and is called the upper invariant 
measure. In other words, let (SiS'(t) denote the probability measure at time t 
for initial probability measure 5i which is the pointmass rj = i(i — 0, 1). Then 
z/ A = lim^oo 5iS(t). Then self-duality of the process implies that v\(A) = 
lim^oo P(£f = 0). The correlation identities for v\(A) can be obtained as 
follows: 



Theorem 1.1 For any A EY , 



E [MAU{y})-u x (A)]+J2[MA\{x})-u x (A) 



A 

x€A cc[=l xGA 



0. 



From now on we consider the one- dimensional case. We introduce the fol- 
lowing notation: 

»x(o) = u x ({0}), u x (oo) = u x ({0, 1}), u x (o x o) = u x ({0, 2}), . . . . 
By Theorem II. 1[ we obtain 
Corollary 1.2 

(1) 2Az/ A (oo) - (2A + l)z/ A (o) + 1 = 0, 

(2) Az/ A (o o o) - (A + l)u x {oo) + u x (o) = 0, 

(3) 2Xu x (o o oo) + u x (o x o) - (2A + 3)z/ A (o o o) + 2u x (oo) = 0, 

(4) Az/ A (o o xo) - (2A + l)u x (o x o) + Az/ A (o o o) + z/ A (o) = 0. 
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The detailed discussion concerning results in this section can be seen in 
Konno [21 El- If we regard A, V\(o), V\(o o o), . . . as variables, then 

the left hand sides of the correlation identities by Theorem 11.11 are polyno- 
mials of degree at most two. In the next section, we give a new procedure 
for getting a series of approximations for extinction probabilities based on 
the Grobner basis by using Corollary 11.21 As for the Grobner basis, see [I] , 
for example. 

2 Our results 

Put X = Z/ A (o), y = Z/ A (oo), Z = V X {o O o), W = V x (o X o), S = V X {° ° 

oo), u = v\(o o xo). Let -< denote the lexicographic order with A -< x -< 
y^w^z^u^s. For m = 1, 2, 3, let I m be the ideals of a polynomial 
ring M[xi, X2, ■ ■ ■ , £n(m)] over E as defined below. Here xi = A, x% = x, X3 = 
y, X4 = z, X5 = w, xq = s, x-j = u and n(l) = 3, n(2) = 4, n(3) = 7. 

2.1 First approximation 

We consider the following ideal based on Corollary 11.21 (1): 

(5) h = ( 2Xy - 2Ax - x + 1, y - x 2 ) C R[A, x, y}. 

Here y—x 2 corresponds to the first (or mean-field) approximation: ^'(00) = 
(z/ A 1} (o)) 2 . Then 

(6) G x = {(x-l)(2Ax-l), y-x 2 } 

is the reduced Grobner basis for I\ with respect to -<. Therefore the solution 
except a trivial one x(= y) = 1 is x = v x (o) = 1/(2A). Remark that the 
trivial solution means that the invariant measure is 5q. From this, we obtain 
the first approximation of the density of the particle, p x = E Vx (r](x)), as 
follows: 



(!) _ 1 „(1) 



2A 



(7) Pr = l-^H°J 2A 

for any A > 1/2. This result gives the first lower bound Ac of the critical 
value A c of the one-dimensional contact process, that is, A^ 1} = 1/2 < A c . 
However it should be noted that the inequality is not proved in our approach. 
The estimated value of A r is about 1.649. 
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2.2 Second approximation 

Consider the following ideal based on Corollary 11.21 (1) and (2): 
7 2 = ( 2Ay — 2\x — x + 1, Xz — Xy — y + x, xz — y 2 ) C R[A, x, y, z). 

Here xz—y 2 corresponds to the second (or pair) approximation: i )^^ (°° 
o) = (^ 2) (oo)) 2 . Then 

G 2 = {(x - 1)((2A - l)x - 1), 1 + 2\(y - x) - x, 

—y — yx + 2x 2 , —z — y(2 + y) + Ax 2 } 

is the reduced Grobner basis for J 2 with respect to -<. Therefore the solution 
except a trivial one x{= y = z) = 1 is x = v\ (o) = 1/(2A — 1). As in a 
similar way of the first approxaimation, we get the second approximation of 
the density of the particle: 

(2) _ 2(A-1) 
Px ~ 2A - 1 ' 

(2) 

for any A > 1. This result implies the second lower bound Ac = 1. We 
should remark that if we take 

I' 2 = ( 2\y — 2Ax — x + 1, \z — \y — y + x, y — x 2 , z — x 3 ) C M[\, x, y, z], 
then we have 

G' 2 = {z-l,y-l,x-l} 

is the reduced Grobner basis for I' 2 with respect to -<. Here y — x 2 and z — x 3 
correspond to an approximation: u^ 2 \oo) = (is^\o)) 2 and ° °) = 

(y^ 2 •'(o)) 3 , respectively. Then we have only trivial solution: x = y = z = 1. 

2.3 Third approximation 

Consider the following ideal based on Corollary 11.21 (l)-(4): 

J 3 = ( 2\y — 2\x — x + 1, \z — \y — y + x, 

2Xs + w - (2A + 3)z + 2y, Xu - (2A + l)w + Xz + x, 

ys — z 2 , xu — yw ) C E[A, x, y, z, w, s, u}. 
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Here ys—z 2 and xu—yw correspond to the third approximation: (oo)u^ (oo 
oo) = (y| 3 '(ooo)) 2 and i/f (o)i^f (oo xo) = (00)1/^ (o x o), respectively. 
Then 

G 3 = {(x - 1)((12A 3 - 5A - l)x 2 - 2A(2A + 3)x - A + 1), . . .} 

is the reduced Grobner basis for I 3 with respect to -<. Therefore the solution 
except a trivial one x = 1 is x = v^\o) = (A(2A + 3) + VD)/(12X s - 5A - 1), 
where D = 16A 4 + 4A 2 + 4A + 1. Then we obtain the third approximation of 
the density of the particle: 

( 8) ,? - 4A ' 3A2 - A - 3 ' 



12A 3 - 2A 2 - 8A - 1 + y/D 

for any A > (1 + \/37)/6. This result corresponds to the third lower bound 
A^ 3) = (1 + V37)/6 « 1.180. 



3 Summary 

We obtain the first, second, and third approximations for the extinction 
probability, the density of the particle, and the lower bound of the one- 
dimensional contact process by using the Grobner basis with respect to a 
suitable term order. These results coincide with results given by the Harris 
lemma (more precisely, the Katori-Konno method, see [3]) or the BFKL 
inequality [5] (see also |3j). As we saw, the generators of I m in Section 2 
have degree at most two in X\,x%, . . ., such as 2\y — 2Xx — x + 1, ys — z 2 in 
the case of I3. We expect that this property will lead to get the higher order 
approximations of the process (and other interacting particle systems having 
a similar property) effectively. 

Acknowledgment. The author thanks Takeshi Kajiwara for valuable dis- 
cussions and comments. 
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